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1.0  Introduction 


In  recent  years,  there  has  been  much  interest  in  exploring  the 
impact  of  learning  and  changes  in  production  rate  on  program  costs. 
Most  researchers  agree  that  learning  is  an  important  determinant  of  cost, 
but  agreement  on  the  cost  impact  of  production  rate  changes  has  been 
less  certain.  Still,  common  sense  and  economic  theory  suggest  that 
production  rate  should  be  an  important  determinant  of  cost.  This 
importance  is  also  suggested  by  the  fact  that  cost  penalties  for  produc¬ 
tion  rate  changes  now  occur  in  some  department  of  defense  contracts  (3). 

This  paper  does  not  present  a  theoretical  justification  for  the 
integration  of  learning  curves  with  traditional  neoclassical  economic 
theory.  The  general  theoretical  framework  for  this  paper  is  published 
in  previous  research  (14,  15).  The  purpose  of  this  paper  is  to  extend 
the  range  of  applicability  of  the  general  framework  by  considering  a 
previously  unexplored  specification.  In  particular,  this  paper  explores 
the  joint  production  situation,  where  learning  and  output  are  simul¬ 
taneously  produced,  and  a  model  is  presented  that  has  potential 
application  in  the  airframe  industry.  The  theoretical  properties  of  the 
model  are  explored,  and  a  cost  minimizing  solution  is  presented. 
Finally,  a  strategy  is  proposed  for  adapting  the  model  to  a  particular 
airframe  program. 


2.0  Background 

The  production  characteristics  of  the  airframe  industry  place 
unusual  demands  on  cost  estimators.  The  situation  is  characterized  by  a 
small  number  of  units  produced,  frequent  design  changes,  and  consider¬ 
able  political  uncertainty.  Both  learning  curves  and  neoclassical  cost 
functions  have  been  used  to  model  this  unique  production  situation,  but 
it  has  only  been  in  recent  years  that  the  two  approaches  have  been 
combined.  Much  of  the  early  work  that  acknolwedged  a  need  to  combine 
the  two  modeling  situations  was  lacking  in  terms  of  theoretical  rigor 
(1,  4,  5).  Also,  many  early  application  oriented  studies  were  empirical 
in  lieu  of  being  firmly  grounded  in  economic  theory  (2,  6).  Most 
recently,  significant  advances  have  been  made  in  modeling  the  made  to 
order  production  situation.  The  first  dynamic  optimization  models  were 
presented  by  Rosen  (7)  and  Washburn  (9),  but  neither  of  these  models  was 
definitive  enough  for  empirical  application.  Recently  a  more  definitive 
model  has  appeared  (14),  and  the  result  has  been  a  series  of  applica¬ 
tions  to  airframe  production  programs  (10,  12). 

The  basic  modeling  framework  was  developed  for  a  firm  producing 
to  an  order  which  specifies  a  quantity  ane  »  "elivery  date  for  output. 
A  neoclassical  production  function  is  augt.’.  with  a  learning  hypoth¬ 
esis,  and  the  cost  of  production  is  minimi*.  to  yield  optimal  time 
paths  of  both  production  and  cost.  The  model  presented  in  this  paper 
adheres  to  the  same  general  theoretical  framework. 
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3.0  The  Model 


Consider  the  situation  where  learning  and  output  are  produced  by 
two  separate  production  technologies.  The  time  path  of  output,  or 
production  rate,  is  presented  as  a  function  of  a  single  variable  com¬ 
posite  resource  and  the  cumulative  stock  of  knowledge.  Also,  the  time 
path  of  experience  rate,  a  decision  variable,  is  presented  as  a  function 
of  the  same  composite  resources  and  the  cumulative  stock  of  knowledge. 
The  assumption  is  that  the  relative  prices  of  the  resources  in  the 
composite  resource  do  not  change,  and  cost  is  measured  in  the  units  of 
the  variable  resource.  The  variables  of  the  model  are  defined  below: 

q(t)  =  output  rate  of  the  program  at  time  t, 

£(t)  =  rate  of  experience  at  time  t. 

t 

Q(t)  =  J  q(T)dt  =  cumulative  output  at  time  t, 
o 


t 

L(t)  =  J  £(t)dx  =  cumulative  stock  of  knowledge  at  time  t, 
o 

Y  =  a  returns  to  scale  parameter, 

^  =  a  returns  to  scale  parameter, 
a  ~  a  learning  parameter, 

6  =  a  learning  parameter, 

C  =  variable  program  cost, 

T  =  time  horizon  for  the  production  program, 

V  =  volume  of  output  to  be  produced  by  T,  l 

M  =  initial  stock  of  knowledge, 

a^  =  constant  term, 

a^  =  constant  term. 

One  possible  specification  is  two  Cobb-Douglas  production  functions, 
that  is, 

q(t)  =  a1xq1/Y(t)L°'(t),  (1) 

and 

£(t)  =  a  x£1/P(t)L6(t). 


By  -- 
D :  " ' 

A 7"  Vi 


l 


(2) 


3 


With  this  specification  the  use  rate  of  the  composite  resource  is 
segregated  into  two  parts,  that  allocated  to  output  x  (t)  and  that 
allocated  to  learning  X«(t).  These  inputs,  combined  with  ?he  cumulative 
stock  of  knowledge  L(tf,  are  used  to  produce  two  products,  output  q(t) 
and  learning  £(t).  The  following  assumptions  define  the  admissible 
ranges  for  the  parameters:  0  <  a  <,  0  <  6  <,  y  >  1,  and  3  >  1. 


The  objective  of  the  firm  is  to  minimize  its  cost  of  production 
subject  to  the  production  function  constraints.  This  may  be  stated  as 
T 


Min  C 


/ 

0 


( x  ( t )  +  x£(t)]dt 


(3) 


subject  to: 

q(t)  =  alXq1/Y(t)La(t),  (4) 
A(t)  =  a2x£1/P(t)L6(t),  (5) 
Q(0)  =  0,  (6) 
Q(T)  =  V,  (7) 
L(0)  =  M,  (8) 


and 


L(T)  *  free. 


(9) 


The  solution  of  (4)  and  (5)  for  x  (t)  and  x^Ct)  yields  the  following 
resource  requirement  functions:  ^ 

Xq(t)  *  qY(t)a^  YL  °Y(t) .  (10) 

and 

x£(t)  =  £P(t)a2'PL‘6p(t).  (11) 

Substituting  the  resource  requirement  functions  into  the  objective 
functional  eliminates  the  production  function  constraints.  The  objec¬ 
tive  functional  is  now  stated  as 
T 

Min  C  =  f  (qY(t)a.’VBY(t)  +  £P(t)a  ’PL’6P(t) jdt.  (12) 

0  l 

A  transformation  of  the  objective  functional  simplifies  the  solution 
procedure.  Let 

Z(t)  =  L1"6(t)/(l-6).  (13) 


This  implies  that 
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z(t)  =  dZ/dt  =  L"°(t)t(t).  (15) 

After  fluking  the  appropriate  substitutions  the  tranfortned  problem  is 

Min  c  =  ;V(t)a 

0  ft  «  (16) 

♦  2P(t)a2‘P)dt 

subject  to: 

Q(0)  =  0,  (17) 

Q(T)  =  V,  (18) 

Z(0)  *  M^/Cl-d),  (19) 


Z(T)  =  free.  (20) 

An  equivalent  way  to  present  the  above  problem  is  as  a  problem  in 
optimal  control  theory.  The  objective  is  stated  as 


Min  C  =  /  (u^(t)a  -V^/(1-6)(l-6)'^/(1-6) 
0  1  i 

*  ..  P  -Pn.. 


+  u2  (t)a2  ™ 


subject  to: 


q(t)  =  u^t), 
z(t)  =  u2(t), 

Q(0)  =  0, 

Q(T)  =  V, 

Z(0)  =  M1"e/(l-«), 


Z(T)  =  free.  (27) 

The  control  variables  for  the  problem,  u^(t)  and  u2(t),  are  the  time 
rates  of  change  of  the  state  variables,  i.e.,  u^(t)  -  q(t)  and  u2(t)  = 
z(t).  The  Hamiltonian  function  for  the  problem  is 
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H  =  u1Y(t)a1_V0,'!f/a"6)(t)(l-6)‘0,Y/{1-6)  + 

ug^(t)a2  ^+X1(t)u1(t)+X2(t)u2(t) .  (28) 

The  necessary  conditions  for  defining  an  extremal  require  that  the 
equations  of  motion,  the  adjoint  conditions,  and  the  Hamiltonian  condi¬ 
tions  hold  simultaneously.  The  equations  of  motion  are 

aH/3^  =  q(t)  =  Ul(t),  (29) 

3H/3X2  =  z(t)  =  u2(t).  (30) 

The  adjoint  conditions  are 

dXx/dt  =  -3H/3Q  =0,  (31) 


dX2/dt  =  -3H/3Z  =  [aY/(l-6)Ju1Y(t)a1'Y 

zM-«v-i)(t)(1.S)-n/«-a).  <32) 

The  Hamiltonian  conditions  are 

3H/3u.  a  Yu1Y‘1(t)a1_YZ-0'Y/a"6)(t) 

(33) 

(l.6)-aY/(l-6)  +  Aj(t)  .  Q> 

3H/3u2  *  pu2P*1(t)a2“P  +  X2(t)  =  0.  (34) 

The  simultaneous  solution  of  these  conditions  requires  solving  two 
second  order  nonlinear  differential  equations.  This  implies  that  there 
are  four  constants  to  be  determined  by  the  boundary  conditions.  Three 
of  the  constants  are  determined  by  the  given  boundary  conditions,  and 
the  fourth  is  given  by  the  natural  boundary  condition,  that  is,  since 
Z(T)  is  free,  the  condition  z(T)=0  determines  the  fourth  constant. 


The  differential  equations  that  follow  from  the  necessary  con¬ 
ditions  are  the  Euler-Lagrange  equations  of  the  calculus  of  variations, 
i.e. , 


M0-l)zS‘2(t)a2'Vz/dt2)  *  lay/d-SjlqV) 

(35) 

and 

(360 

where  kj  is  a  constant  of  integration. 

The  simultaneous  solution  of  the  Euler-Lagrange  equations  yields 
the  optimal  time  paths  of  both  production  and  experience  rates.  The 
solution  proceeds  as  follows.  Solve  (36)  for  q(t)  and  state  the  result 
in  compact  notation  as 
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q(t)  =  k1Z°(t)  (37) 

where  rpay/ ( (1-6)  (y-1)  ]  ■  This  expression  for  q(t)  may  be  substituted 
into  (35),  and  the  number  of  equations  is  reduced  by  one.  The  single 
necessary  condition  is  stated  as 

d2Z/dt2  =  p'10-l)'1z2'P(t)a2P{aY/(l*6)J 

z(l+ay+6Y-Y-6)/((Y-l)(l-6)](t)kiY(1.6)nai*Y  (38) 

This  is  a  second  order  nonlinear  differential  equation  which  may  be 
stated  in  compact  notation  as 

d2Z/dt2  =  Az2"P(t)Zn'1(t)  (39) 

where  the  constant  term  is 

A  =  -p'1(p-l)'1a2PlaY/(l-6)U1Y(l-6)"0'Y/(1-6)a1"Y  (40 


and 


n-1  =  (l+QfY+<5Y-Y“<5)/Ky-l)(l-6)J. 


(41) 


After  the  series  of  transformations  to  achieve  a  reduction  in  order,  the 
following  differential  equation  is  obtained: 

dZ/dt  =  [APZn(t)/q  +  k2J1/P.  (42) 


This  function  may  be  inverted,  leaving  t  as  the  following  function  of 
Z(t): 


.  .  JZ(t)[A0Z%  +  kj'1/PdZ  +  k 

t  =  1 


(A3) 


An  expression  is  needed  that  gives  Z  as  a  function  of  t.  This  allows 
you  to  apply  the  boundary  conditions  on  Z(t),  determine  the  constants, 
and  define  one  of  the  extremals.  Since  (37)  relates  Z(t)  to  q(t),  you 
should  be  able  to  apply  the  boundary  conditions  on  Q(t)  and  define  the 
second  extremal.  Unfortunately  this  procedure  is  easily  stated  but  not 
easily  executed. 


The  next  step  in  the  solution  procedure  is  just  a  restatement  of 
(43)  in  a  slightly  different  form.  After  some  algebraic  manipulation 
(43)  may  be  written  as 


=  k2‘l/PM1_6/(l-6) 


/Z(t)  [  1  +-A0Zn/rtk, f 1 7  PdZ+k, . 


(44) 


The  value  of  k,  is  determined  by  returning  to  (40) ,  which  allows  you  to 
apply  the  natural  boundary  condition  z(t)  -  0.  The  resulting  expression 
for  k2  is 

k2  =  -A0Zn(T)/n, 


and  (44)  may  be  written  as 


(45) 


.  .  [-A3Zn(T)/n)"1/P  JZ(t)[l-z‘n(T)Znj'1/P  dZ+k_ .  (46) 

t  =  M1  °/(l-6)  3 


This  integral  appears  uninviting,  but  a  change  of  variables  leads  to  a 
solution.  Let 

y  =  Z"n(T)Zn.  (47) 


The  integral  may  now  be  restated  as 

t  =  i-A3zn(T)/nr1/pz(T)n"1 

Z'n(T)Zn(t) 

Z"n(T)(M1'5/(l-6)]n 

,,  ,u-i/p)-i  i/n-i,  „  . 

( 1-y)  r  y  1  dy  +  k3 

which  is  a  form  of  the  incomplete  beta  function. 


(48) 


The  solution  is  complete  except  for  the  determination  of  the 

constants  of  integration.  At  this  point  the  value  of  k.  is  unknown,  and 
there  is  also  an  unknown  constant,  k.  ,  in  the  A  term.  To  determine 
these  constants,  first  notice  that  (26;  implies  that  1^=0.  The  deter¬ 
mination  of  the  second  constant  requires  a  little  more  effort.  The 

strategy  is  to  find  an  expression  for  [-ApZr,(T)/r|]  that  is  in  terms  of 
known  constants.  The  solution  is  as  follows.  Equation  (42)  may  be 
written  as 

dt  =  [-ApZ‘n(T)/nf 1/P  ll-Z"n(T)Zn(t)]'1/PdZ.  (49) 

It  also  follows  from  the  Euler-Lagrange  equations  that 

Q(T)  *  k  JtZ(x)dt.  (50) 

0 


After  evaluating  (49)  at  t  and  changing  variables  in  (50),  it  is  pos¬ 
sible  to  write  Q  as  a  function  of  Z,  i.e., 

q(z)  *  k1(-Apzn(T)/n]"1/p 

Z(t)  „  „  n  I/O 

.  *  J  Zn(l-z’n(T)Zn]  /PdZ.  (51) 

H^/tW) 

This  integral  results  in  an  expression  that  is  suitable  for  applying  the 
boundary  conditions  on  Q(t).  In  other  words,  (50)  is  transformed  into 
an  expression  that  is  integrated  with  respect  to  Z  in  lieu  of  T.  Con¬ 
tinuing  with  the  solution,  let 

R  =  k1[-ApZn(T)/nl"1/P,  (52) 


and 


y  =  Z-n(T)Zn. 


(53) 
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Equation  (51)  may  now  be  stated  as  a  form  of  the  incomplete  Beta 
function.  The  appropriate  integral  is 


Q(y)  =  RZn+1(T)n_1 
Z~n(T)Zn(t) 


Z*n(T) [M1’6/ (1-6)] 


/  y^/n+l)-l(1.y)(l-l/P)-ldy  +  k 


(54) 


It  is  now  possible  to  apply  the  boundary  conditions  on  Q.  The  initial 
condition,  Q(0),  =  0,  implies  that  k4  =  0.  The  final  condition,  Q(T)  = 
V,  implies  a 


R  =  VZ‘1_n(T)n  {  S  y(1/n+1)"1(l-y)(1‘1/0)'1dy}'1 


Zn(T) (M1-6/ (1-6) ]n 


which  is  another  form  of  the  incomplete  Beta  function. 


Since  all  of  the  integration  constants  are  known,  it  is  possible 
to  write  an  expression  that  links  optimal  Z  with  T.  The  expression  is 

z“n(T)Zn(t) 

t  =  ^  jy1/n-x{1.y)(1-1/P)-idy  (56) 

Z"n(T) [M1_6/(l-6)]n 


where  c.  is  comprised  completely  of  constants.  By  using  the  inverse 
incomplete  Beta  function  in  (56),  it  is  possible  to  determine  Z(t)  for 
any  t.  Optimal  Z(t)  determines  optimal  L(t)  via  (13),  and  optimal  Z(t) 
determines  optimal  Q(t)  via  (37).  The  theoretical  solution  is  now 
complete,  but  unfortunately  the  model  is  not  in  a  form  that  is  suitable 
for  application  to  most  airframe  programs. 


4.0  Strategy  for  Application 

The  theoretical  model  requires  some  adjustment  before  it  is  in 
the  proper  form  for  application  on  most  airframe  programs.  The  problem 
is  that  the  model  is  not  stated  in  terras  of  variables  that  are  observ¬ 
able  or  even  measurable.  There  is  no  convenient  way  to  measure  the 
current  stock  of  knowledge,  and  there  is  even  some  doubt  about  the 
proper  way  to  measure  production  rate  (8).  However,  there  are  two 
quantities  that  are  usually  reported  with  a  degree  of  regularity. 
Direct  labor  hours,  either  by  airframe  or  unit  time,  is  usually  avail¬ 
able,  and  there  is  usually  some  information  on  delivery  schedules  and 
perhaps  lot  sizes.  The  latter  information  is  useful  for  assigning  cost 
to  time  periods  to  develop  a  data  series  suitable  to  estimate  the  model 
parameters.  The  following  procedure  provides  a  method  for  transforming 
the  model  so  that  it  may  be  used  in  applications  where  the  only  avail* 
able  data  is  cost  per  unit  time. 
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The  objective  to  rewrite  the  optimal  inverse  function  (48)  so 
that  t  is  a  function  of  x,  a  quantity  that  is  observable  in  the  data. 
The  total  resource  requirement  function  is  the  sum  of  the  individual 
resource  requirement  functions,  i.e., 

x(t)  =  xq(t)  +  x£(t)  (57) 

After  substituting  (10)  and  (11)  into  (57)  and  using  (15),  the  combined 
resource  requirement  function  may  be  rewritten  as 

x ( t )  =  qY(t)a1“Y2‘0,Y/(1_6)(t)  +  ZP(t)a2‘0.  (58) 

The  strategy  is  to  eliminate  q(t)  and  z(t)  from  the  above  expression. 
This  leaves  an  expression  which  may  be  solved  for  Z(t)  as  a  function  of 
x (t ) . 


The  following  procedure  is  used  to  eliminate  z(t).  Equation 
(42)  implies  that 

zP(t)  =  A0Zn(t)/n  +  kr  (59) 

If  this  result  is  substituted  into  the  resource  requirement  function, 
z(t)  is  eliminated.  The  resources  required  may  now  be  written  as 

x(t)  =  qY(t)a1'YZ"0fY/(1"6)(t)  +  a2‘PApZn(t)/n  +  k^"*.  (60) 

To  eliminate  q(t),  use  the  Euler*Lagrange  equation  (37).  Solve  the 
Euler-Lagrange  equation  for  q(t),  and  substitute  into  (60)  to  obtain  the 
desired  result.  After  solving  for  Z(t),  the  optimal  expression  is 

Z(t)  =  (D+Ef  1/fl[x(t)  -  k2a2‘P]1/n.  (61) 

where  D  and  E  are  functions  of  known  constants.  If  (61)  is  substituted 
as  the  upper  limit  of  integration  in  the  optimal  inverse  function  (48), 
the  result  is  an  expression  that  gives  t  as  a  function  of  x(t)  at  any 
point  in  time,  i.e., 

t  =  z(-o+0-y)/  (0-k)  (T)n(i-pn)  /  (P-y) 

fl*Y/(0*Y)  (fj-j)*/ (P"Y)a  ”P/  (P"V)a  Y/  0"Y) 

(62) 

Z'n(T)(D+E)_1[x(t)-k2a2'P] 

z“n(T)lMl‘a/(l-6)Jrrr 

(1-y)(1-1/P)‘1y1/n_1dy. 


Although  the  transformation  is  complete,  the  function  is  still  not 
appropriate  for  estimation  since  the  quantity  that  is  observable  is  not 

x(t),  but  cumulative  x(t)  over  some  interval,  i.e.,  C  =  /  x(t)dt.  As 


a  final  adjustment,  (61)  is  solved  for  x(t)  and  integrated 
interval  of  time.  The  required  integration  is 


over  an 
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C=X(t1)-X(t0)  =  S  1[(D+E)Zn(t)  +  k2a"P]dt 


(63) 


If  the  Euler-Lagrange  equation  (37),  is  used  to  eliminate  Z(t)  from 
(63),  the  integral  becomes 


C  = 


;  1[(D+E)k  _1q(t)  *  k  a  -P]dt. 
*0 


(64) 


After  performing  the  integration  and  applying  the  boundary  condition  on 
Q(t),  the  final  result  is 

C=  (D+E)k1'1(Q(t1)-Q(t0)]+k2a2'P(t1-t0).  (65) 

This  result  is  the  basic  estimable  relationship  of  this  paper.  Since 
Q(t)  is  known  for  any  value  of  Z(t)  by  equation  (51),  and  Z(t)  is  known 
for  any  value  of  t  by  equation  (48),  X(t)  is  known  for  any  value  of  t. 


Previous  experience  with  estimating  parameters  in  this  type  of 
model  suggests  that  in  many  cases,  the  model  is  over  parameterized  (12). 
Our  future  papers  will  explore  the  economic  implications  of  a  simpler 
reparameterized  version  of  (65).  The  specification  is 

C  =  C^QCc^-QUq)]  +  c2(t1-tQ)  (66) 

The  complex  expressions  for  Q(t.)  provide  the  necessary  functions  for 
estimating  the  model  parameters. 


Simulations  with  (66)  indicate  that  this  model  generates  total 
cost  curves  whose  shape  are  very  similar  to  those  encountered  in  prac¬ 
tice.  Figure  1  shows  various  simulated  cost  curves  using  apriori  para¬ 
meter  estimates  and  the  volume /delivery  schedule  values  (V  and  T)  from 
the  F102  airframe  program.  Notice  that  the  model  is  very  responsive  to 
changes  in  the  returns  to  knowledge  parameter  p.  Also,  Figure  2  gives 
simulations  of  the  cumulative  stock  of  knowledge  using  these  same  values 
of  p.  This  curve  is  extremely  hypothetical  since  arbitrary  values  were 
assigned  to  the  terminal  stock  of  knowledge.  Our  current  research 
includes  the  estimation  of  the  parameters  in  this  model  using  nonlinear 
least  squares. 


5.0  Conclusion 

This  paper  examines  a  theoretical  model  of  the  "made-to-order" 
production  situation.  The  model  is  unique  in  that  it  examines  the  case 
where  output  and  knowledge  are  produced  by  two  separate  production 
technologies.  The  solution  for  the  model  yields  an  optimal  estimable 
relationships  which  has  potential  for  application  in  the  airframe 
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industry.  The  usefulness  of  the  model  will  be  tested  in  a  later  paper 
where  the  model  is  applied  to  the  F102  airframe  program. 
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